Long Range Coulomb Interactions and Nanoscale Electronic Inhomogeneities in 

Correlated Oxides 



Vijay B. ShenoyiE Tribikram Gupta^Il] H. R. Krishnamurthyi^2|i] g^^^j V. Ramakrishnan^.i^^ 
^Centre For Condensed Matter Theory, Indian Institute of Science, Bangalore 560 012, India 
^ Jawaharlal Nehru Centre for Advanced Scientific Research, Jakkur, Bangalore 560 064, India 
"^Department of Physics, Banaras Hindu University, Varanasi 221 005, UP, India 

(Dated: April 12, 2009) 

Electronic, magnetic or structural inhomogeneities ranging in size from nanoscopic to mesoscopic 
scales seem endemic, and are possibly generic, to colossal magnetoresistance manganites and other 
transition metal oxides. They are hence of great current interest and understanding them is of 
fundamental importance. We show here that an extension, to include long range Coulomb interac- 
tions, of a quantum two-fluid £ — b model proposed recently for manganites [Phys. Rev. Lett., 92, 
157203 (2004)] leads to an excellent description of such inhomogeneities. In the £ — b model two 
very different kinds of electronic states, one localized and polaronic {£), and the other extended or 
broad band (fe) co-exist. For model parameters appropriate to manganites, and even within a sim- 
ple dynamical mean-filed theory (DMFT) framework, it describes many of the unusual phenomena 
seen in manganites, including colossal magnetoresistance (CMR), qualitatively and quantitatively. 
However, in the absence of long ranged Coulomb interaction, a system described by such a model 
would actually phase separate, into macroscopic regions of I and b electrons respectively. As we 
show in this paper, in the presence of Coulomb interactions, the macroscopic phase separation gets 
suppressed, and instead nanometer scale regions of polarons interspersed with band electron puddles 
appear, constituting a new kind of quantum Coulomb glass. We characterize the size scales and 
distribution of the inhomogeneity using computer simulations. For realistic values of the long range 
Coulomb interaction parameter Vb, our results for the thresholds for occupancy of the 6 states are in 
agreement with, and hence support, the earlier approach mentioned above based on a configuration 
averaged DMFT treatment which neglects Vb; but the present work has new features that can not 
be addressed in the DMFT framework. Our work points to an interplay of strong correlations, 
long range Coulomb interaction and dopant ion disorder, all inevitably present in transition metal 
oxides, as the origin of nanoscale inhomogeneities rather than disorder frustrated phase competition 
as is generally believed. As regards manganites, it argues against explanations for CMR based on 
disorder frustrated phase separation and for an intrinsic origin of CMR. Based on this, we argue 
that the observed micrometer(meso)-scale inhomogeneities owe their existence to extrinsic causes, 
eg. strain due to cracks and defects. We suggest possible experiments to validate our speculation. 

PACS numbers: 75.47.Lx, 71.30.-f-h, 71.38.-k, 75.47.Gk 



I. INTRODUCTION 



In the last decade or so, a number of experiments on 
several families of transition metal oxides have provided 
evidenceii^ that these often consist of patches of two 
different kinds of electronic/structm-al/magnctic states. 
The patches range in size from nanometers to microm- 
eters, and can be static or dynamic. In doped man- 
ganites (Re(i_a;) AkjjMnOs virith Re and Ak being rare- 
earth and alkaline earth ions respectively) where this 
phenomenon seems most widespread, one can have in- 
sulating, locally lattice distorted (in some cases charge 
ordered) regions coexisting with metallic, lattice undis- 
torted one s i^iii^i^ili^ i ^i^^^ii cuprates, an antiferromag- 
netic insulating state and a metallic (or superconducting) 
state seem close in energy and may coexist under some 
conditions, eg. as nanoscale stripes (see, for example, the 
references ini^). There is a view that the proximity of 
two states with very different long range orders (LRO) 
is a defining characteristic of these systems. For exam- 
ple, such a view has been most forcefully put forth in the 
work of Dagotto and coUaborators^ii^iiiii^ on the basis 



of numerical simulations of finite sized samples of simpli- 
fied lattice models capturing one or more key features of 
manganites. In Ref.lla. models with Mn Cg electrons hav- 
ing Jahn- Teller coupling to classical phonons and double 
exchange coupling to Mn t2g core spins and antiferromag- 
netic interactions Jaf between the core spins are studied 
for x=0.25 on lattices of size up to 12 x 12. The clean 
system has a first order transition (ending at a bi-critical 
temperature from a ferromagnetic metallic (F-M) phase 
for small Jaf to an antifcrromagnctic(AF-I). charge or- 
dered, insulating phase when Jaf crosses a threshhold 
JAF-t- The magneto- resistance exhibits CMR like fea- 
tures when Jaf is tuned to be in the vicinity of the 
transition coupling Jaf-i- But the tuning becomes less 
of a requirement in the presence of disorder, which fur- 
ther generates nanoscopic interpenetrating patches of the 
two regions. These arc then argued to be generic fea- 
tures responsible for (nanoscale) two 'phase' coexistence, 
colossal magnetoresistance (CMR) and other phenom- 
ena observed in manganites. Larger length scale (mi- 
cron scale) coexistence is also attributed to disorder ef- 
fects amplified by proximity to a critical value of compel- 
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ing parameters However, the existence of manganites 
with neghgible frozen disorder but showing CMR and 
other characteristic phenomena— suggests that these ef- 
fects are intrinsic, and that nano- and micron-scale inho- 
mogeneities could arise from other, intrinsic, causes. 

In particular, the Coulomb interaction which is in- 
evitably present can play a significant role if the two 
phases involved have different charge densities, and will 
suppress phase separation, leading to nanoscopic elec- 
tronic inhomogeneity or 'phase' coexistence. This is the 
ubiquitous effect we investigate here quantitatively, for 
the first time, in conjunction with the two fluid £ — b 
mode l^^i^^'^^ proposed recently for manganites. This 
model invokes two very different types of electronic de- 
grees of freedom, one polaronic and localized (called £) 
with associated Jahn- Teller lattice distortion, and the 
other (called b) forming a broad band, with no lattice 
distortion, and moving around primarily on sites not oc- 
cupied by the polarons. It has been show n^^'^^ that this 
approach explains much of the unusual behavior of man- 
ganites qualitatively and quantitatively. The work here 
shows in detail how the inclusion of electrostatic Coulomb 
interactions in the £ — b model leads to nanometer scale 
inhomogencities, and quantifies this effect. More gener- 
ally, our results describe the intrinsic emergence of such 
nanoscale electronic inhomogencities in any many elec- 
tron system with two very different microscopic states of 
comparable energy. The model is different from disorder 
based phase separation or domain formation as seen in 
computer simulations of simple models^ or in statistical 
Imry-Ma— type arguments. The frequently observed 
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micron scale inhomogeneity in many manganite; 
probably related to strain effects, which are long ranged 
and unscreened. There are no calculations of this ef- 
fect, though simulations of elastic strain effects and inho- 
mogencities in Jahn- Teller distorted lattice systems have 
been made^l on simplified models. 

Specifically, in the ^—6 model, at each lattice site there 
can be two types of electrons, namely a nearly localized 
polaronic one (called £) with site energy —Ejt and a 
broad band one (called 6, with site energy zero (0) and 
hopping amplitude t). There is a strong local repulsion 
U (we consider U ^ oo in this work) present between 
the £ polaron and the b electron. The physical origin and 
parameters of this model for manganites are mentioned 
below (Section III AD and have been described in refer- 
ences Il7lll8lll9ll22l |. If long range Coulomb interactions 



are neglected, the system described by the model phase 
separates, i. e., the ground state consists of two separate 
phases, one a macroscopic region entirely of £ polarons 
at every site, and the other consisting of occupied b band 
states, the b density being fixed by the requirement of 
uniformity of the chemical potential /i, whence ^ must 
equal —Ejt- However, the two phases have very differ- 
ent charge densities. With respect to ReMnOs, which has 
one Cg electron per site, as the reference state, the pola- 
ronic phase is neutral, while the b phase has a positive 
charge density of 1 — fib per site. Furthermore, the doped 



alkaline earth ions (most likely randomly distributed on 
the lattice) supply unit negative charge at each alkaline 
earth site. The relevant "extended £b model" which in- 
cludes long range Coulomb interaction involving these 
charges is described in Section fllBI It leads to the sup- 
pression of macroscopic phase separation and the forma- 
tion of nanoscopic "puddles" of b electrons surrounded by 
regions with only I polarons. We first present a simple 
analytical estimate of the size of these regions, in partic- 
ular their dependence on Vq, the energy parameter that 
determines the strength of the long range Coulomb in- 
teraction (Section III Cp . Next, we describe simulations 
that we have carried out on finite systems of size up 
to 20 X 20 X 20, with the energy levels of the b elec- 
tron puddles calculated quantum mechanically, and the 
Coulomb effects treated in the Hartree approximation. 
The method for the ground state determination is de- 
tailed in Section fll Dl The results of the simulations are 
described and discussed in Section IIIII The £ polarons 
are shown to form a Coulomb glass^^. Electronic states 
in the b electron regions are occupied (up to the chemi- 
cal potential /x) for a critical hole concentration x > Xd] 
these puddles connect percolatively and the system is a 
metal for x > Xc2 > Xci ■ We analyze and discuss in 
detail Xd, Xc2 as well as b clump size distributions as a 
function of Vq , Ejt and the arrangement of the alkaline 
earth ions. We exhibit several examples of real space in- 
homogeneous structures. The spatial autocorrelation of 
holes (absence of polarons), and of holes with alkaline 
earth ions is also elaborated. We compare the effective 
b bandwidth obtained from simulations which have the 
intrinsic nanoscale inhomogeneity mentioned above with 
results from sing le site DMFTii neglecting long range 
Coulomb interactions but forcing homogeneity ('annealed 
disorder'), and find very good agreement. 

In the final part of this paper (Section HV)) we discuss 
some of the approximations made, eg. the assumption of 
random distribution of dopant ions, the effect of disorder 
and its modeling. We also mention a number of impli- 
cations of our results for manganites such as Coulomb 
glass behavior and other signatures of nanoscale inhomo- 
gencities on transport properties. We place our results 
in the context of the inhomogencities observed in experi- 
ments, and suggest some new experiments to investigate 
long range strain effects. A short description of this work 
has been publishedi^ 



II. MODEL AND METHOD 

This section contains four parts. In the first part 
(Sec. ESI, the £b mode]i2J^ is briefly described. The 
second part (Sec. IIIB[) contains a description of the ex- 
tended £b Hamiltonian. A simplified analytical solution 
for the ground state of the extended £b Hamiltonian is 
presented in the third part (Section III Cp . Finally, the 
determination of the ground state of the extended £b 
Hamiltonian is described in the final part (Sec. Ill D|) . 
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FIG. 1: (color online) "Macroscopic phase separation" in the 
£b model in absence of long ranged Coulomb interaction. The 
lighter regions correspond to holes that form a large clump 
(with occupied b states), the darker regions are occupied by 
£ polarons. This simulation is performed with a cube of size 
16x16x16. 



A. Summary of the lb model 

The Mn^"*"'^"*" ions in Re(i_2,)Aka;Mn03 form a simple 
cubic lattice (with lattice parameter taken as a, lattice 
sites indexed by i). Each manganese ion experiences the 
octahedral environment of oxygen atoms and thus the 



d states are crystal field split into t2g and Cg orbitals. 
For the purposes of describing the low energy physics of 
manganites one can treat the t2g levels as always occu- 
pied by three electrons at every site, with parallel spin 
because of a strong Hund's rule, and hence replace them 
by a S" = I , - core spin" , and in the rest of this 
paper we approximate this as a classical spin, SCl where 
51 is a unit vector. The remaining (1 — x) electrons per 
site move around among the eg orbitals with an average 
hopping amplitude t. However, they have a strong Hund 
coupling Jff with the core spin, and there is a strong 
on-site Coulomb repulsion U between the eg electrons. 
Furthermore, they have a strong Jahn- Teller interaction 
with the modes of distortion of the oxygen octahedron 
surrounding the Mn'^"'' ions. As a consequence of this 
and phonon dynamics, two types of low energy, effec- 
tive, eg electronic states called £ and b emerge. The 
polaronic £ state is associated with a local lattice dis- 
tortion in which an eg electron can get self trapped, with 
a binding energy Ejt- The hopping of the ^-polaron 
gets suppressed exponentially by the Huang-Rhys factor 
rj (= (,-EjT/2hujo ^ 1/200, where luq is the frequency of 
the local lattice distortion). The reduction arises from 
the fact that the hop of the polaron from one site to a 
neighboring one, in addition to the transfer of the elec- 
tron, involves the relaxation of the lattice distortion at 
the original site and effecting a similar lattice distortion 
at the neighboring site, and the exponentially small over- 
lap of the corresponding phonon wavefunctions. On the 
other hand, as t >> wq, Cg electrons can also hop fast 
among empty, undistorted or weakly distorted sites with 
essentially the bare amplitude t, leading to the b states. 
The greatly diminished hopping of the polarons allows 
for a useful approximation: at temperatures much larger 
than rit, the polarons may be considered as static. The 
effective £b Hamiltonian can therefore be written as 



i.(T <ij> « 

+U^niinib - fi^{nii + na) - {Jvde - Jse) ^ Si ■ Sj. 

i i 

I 



(1) 



Here sa^nu) and Sihinn,) are spin (number) operators 
corresponding to the £ and b degrees of freedom at site 
i, and fi is the chemical potential determined from the 
condition that (n/ -I- nb) = 1 — x. The last term con- 
tains the "virtual double exchange" ferromagnetic cou- 
pling with strength JyDE^x{l — x)^j^. This coupling 
between the neighboring core spins arises from fast vir- 
tual hopping processes of the ^-electron to a neighboring 
vacant site. Jse is the antiferromagnetic superexchange. 
The above Hamiltonian treats all the main energy scales 



that govern manganite physics. The three largest local 
(strong correlation) energies are the on-site Coulomb or 
Mott-Hubbard repulsion C/ sa 5 eV, the Hund's rule fer- 
romagnetic exchange coupling Jh ~ 2eV and the Jahn- 
Teller energy Ejt (~ 0.6 to 1 eV). The other important 
interactions are the nearest neighboring hopping t sa 0.2 
eV, the doping dependent virtual double exchange, and 
the superexchange Jse ^ 0.01 eV. 

Several simplifying approximations have been made in 
writing the Hamiltonian ([T|). The first is that the ki- 
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netic energy in the Hamiltonian is "orbitally averaged", 
i. e., the hopping ampUtude t represents an average over 
the possible orbital configurations at the two pertinent 
sites. Second, there are no cooperative/long range lat- 
tice effects, i. e., no intersite polaron correlations. Third, 
the virtual double exchange has been approximated as a 
homogeneous interaction Jvde ; in reality it acts only 
between pairs of sites where one is occupied by a polaron 
and the other is empty. Fourth, the Jahn- Teller interac- 
tion has been included only to the extent that it leads to 
the formation of polarons. But its merit is that all the 
important energy scales governing manganite physics are 
included in ([T]), at least in an approximate fashion. Much 
of the previous work on manganites is based on simpli- 
fied models^ '^^i^^ that neglect one or more of these energy 
scales. 

The Hamiltonian ([T]) closely resembles the Falicov- 
Kimball model (FKM)21; in fact, for Jh — »■ oo, and at 
T = 0, when all the spins are completely ferromagneti- 
cally aligned, it is the same as the FKM. The model was 
solvediiii^ii^ using the dynamical mean field theory^^. It 
is successful in capturing the colossal magnctoresistance 
effect, the ferro- magnetic insulating state found in low 
band width manganites, the systematics of the the role 
of R and Ak ion radii, etc. The key point is that for small 
doping X the majority of the Cg electrons get localized as 
i polarons. Because of the strong on-site Coulomb in- 
teraction U (the largest energy scale in the problem) the 
h electrons avoid the polaronic sites, and this reduces 
the effective half bandwidth of the h states from its bare 
value Dq to i'e// which is strongly (x) dependent. At 
zero temperature with f7, Jh oo, and a ferromagnetic 
ordering of the core spins, the effective half bandwidth is 
given hy^ 

Deff^V^Dn. (2) 

Thus, at low doping, the effective band bottom is above 
the polaronic energy level —Ejt and therefore the h band 
is unoccupied and the chemical potential /i is pinned at 
—Ejt- Clearly, beyond a critical doping Xc given by 

the 6-states begin to be occupied, leading to a insulator 
to metal transition. Furthermore, for x > Xc, starting 
from the metallic state at T = 0, the effective band- 
width reduces from its zero temperature value with in- 
crease of temperature. This is because the hopping of 



the h electrons is strongly inhibited on account of their 
Hund coupling Jh to the thermally disordered core spins. 
Hence one gets a thermally induced ferro-metal to para- 
insulator transition when the h band moves above Ejt- 
For T near Tc, the application of an external magnetic 
field causes the core spins to align and thus increases the 
effective bandwidth, leading to increase in the number of 
thermally excited carriers by orders of magnitude. This 
causes colossal magnctoresistance. 

The work cited above did not address the issue of elec- 
tronic inhomogeneities. In fact, an effective homogeneous 
state was assumed in the dynamical mean-field solution. 
This is a drastic assumption, because from previous work 
on the Falicov-Kimball model^S (and as confirmed by our 
computer simulations (see FIG.[T])the ground state of the 

Hamiltonian ([1]) (for the parameter ranges discussed 
above) is known to be a macroscopically phase separated 
state due to the strong on-site Coulomb correlation U- 
AU the polarons cluster on one side of the box; this al- 
lows band states (6) to optimize their kinetic energy by 
moving among the vacant 'hole' sites in the other part of 
the box. But the two portions have drastically different 
electron densities. The motivation for the homogeneity 
assumption made in the DMFT work was that the phase 
separation will of course be prevented by the long ranged 
Coulomb interactions which are always present in the real 
system. In this work, we extend the £b model to include 
long ranged Coulomb interactions to address the issue of 
electronic inhomogeneities, and to check to what extent 
the DMFT homogeneity assumption is valid. 



B. The extended ib Hamiltonian 

The Hamiltonian is most conveniently developed in 
terms of 'hole' operators hi = £i, which create vacant 
sites not occupied by polarons. A hole carries a positive 
unit charge. The charge is counterbalanced by the Ak"*"^ 
ions which have a negative unit charge. The Ak ions 
occupy the sites (^, ^, ^) where l,m,,n are odd integers 
(mimicking the perovskitc type structure of manganites); 
the number of Ak ions is equal to xN {N is the number of 
sites in the model) and they are placed randomly on the 
lattice sites indicated. In the orbital liquid regim e^^'^^i^^ 
with large Hund coupling and on site Coulomb repulsion 
discussed above, and at T = 0, corresponding to a fully 
ferromagnetic alignment of the core spins (spin drops out 
of the problem in this limit), the extended £b Hamiltonian 
for fixed particle number is given by 



ntf = i^Ak + E '^^'i^ + ^" E — ^^'^^ - * E (^l^j- + ^h) + E /^l^^ (4) 



I 



where is the charge operator at hole site i (in units of |e|). i?Ak 

- _ , ti, ;,t!, is the Ak-Ak electrostatic interaction energy (a fixed 
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number for a given realization of Ak distribution in the 
lattice), is the Coulomb potential at the hole site i 
due to Ak ions, Vq is the strength of the Coulomb in- 
teraction between nearest neighbor holes (of the order of 
0.02eV, Vq « O.Olt — O.lt, see discussion below), is 
the distance between holes i and j in units of the lat- 
tice constant. Each hole has an energy penalty of Ejt- 
The sum is over all pairs of hole sites, while (ij) 

denotes a sum only over nearest neighbor hole sites. The 
site charge operators satisfy the constraint 



xN. 



(6) 



The U ^ oo constarint implies that the b electrons can 
only move among 'hole' sites, where hlhi = 1. 

The Hamiltonian (|4]) embodies the competition be- 
tween long range Coulomb interaction which works to 
keep the holes as far apart as possible, while the ki- 
netic energy of band electrons promotes the formation 
of "clumps" . A clump is a collection of hole sites, each 
member of which can be reached from any other mem- 
ber via a sequence of nearest neighbor hops which visit 
only the members of the clump. Thus every configura- 
tion of holes can be broken up into a set of clumps where 
each hole belongs to only one clump. Site delocalized 
electronic b states are possible in a clump that has more 
than one member (a one member clump is called a "sin- 
gleton" ) . The kinetic energy gain can possibly promote 
an £ polaron to occupy the b state in the clump lead- 
ing to an electron "puddle" , which creates an additional 
hole in the system at the site from which the £ polaron 
is removed, which changes the electrostatic energies, etc. 
Thus the number of holes and number of band electrons 
is not individually conserved, but only the constraint © 
is satisfied. The key question is: what is the ground 
state of the Hamiltonian for a realization of the random 
distribution of the Ak ions? The determination of this 
ground state requires an optimization of the number of 
holes and their distribution [or equivalently, the configu- 
ration of the on site £ polarons, which exist at sites where 
the holes are not there], as well as the clump structure, 
the associated b states and their occupancy (constitut- 
ing the 'b' puddles surrounded by the £ plarons), so as 
to achieve the lowest energy of the Hamiltonian An 
exact solution of this problem is beyond the available 
techniques of correlated electron theory, so approxima- 
tions have to be resorted to. In this paper we explore 
two alternate ways; the first one is essentially analytic 
(and very approximate), while the second one is a full 
scale numerical treatment. 



C. A simplified analytic treatment 

A simple, analytical approximation for the ground 
state becomes possible if we assume that the Ak ions 
are distributed homogeneously in space, i. e., model the 
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FIG. 2: Schematic of the ground state used in the analytical 
calculation of Section III CI "Phase separation" takes the sys- 
tem to two distinct type of regions called the "bulk" (regions 
with I polarons) and the clump. The dashed lines indicate the 
assumed periodic nature of the clump distribution - clumps of 
size R (volume R^) are assumed to be arranged in a periodic 
fashion (period L) with intervening "bulk" regions. Site delo- 
calized electronic states are found in the clumps. The initial 
hole density is x (also equal to the background negative charge 
density). The charge density in the clump is [1 — Uh — x) and 
that in the "bulk" is rih — x. nt is the fraction of electrons 
that are promoted to delocalized states, and Uh is the fraction 
of holes created in the 'bulk'. The fractions x, rib and rih are 
related via charge balance. 



charge of the Ak ions as a "jellium" . The approxima- 
tion involved is a variational calculation with an assumed 
ground state as shown in FIG.[21 where clumps of size R 
(volume R^) are taken to be spaced periodically with a 
spacing L. The regions between the clumps - the "bulk" , 
is taken to have concentration of holes rih (the concen- 
tration of holes in the clump is, of course, unity). Charge 
conservation gives 



(1 - nb)R^ + nh{L^ - R^) = xL^ 



whence 



nh 



a^x+ {b- 1) 
a3 - 1 



(7) 



(8) 



with a = L/R. Clearly, the charge density pc in the 
clump is (1 — rih — x) and that in the bulk pB is rih — x, 
I. e., 



PC = 1 ~ - X, pb ^ Uh - X ^ -(1 - Ub - x) 



1 



a3 - 1 



.(9) 
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Associated with this charge distribution, there is an 
electrostatic energy (per "unit ceU" of volume Li^ = 
a^i?^) given by 



The polaron energy (again for a spherical clump) is 



Essinb, a, R) = Vo£{nb, a, R/a) 



(10) 



47r 

Ep{nb,a,L) = —EjrUbR ■ 



(16) 



where a is the lattice parameter of the underlying atomic 
lattice, £{nb,a) is the "Madelung function" accounting 
for the net electrostatic energy in the system. Since a 
total number of UbR^ polarons have been promoted to 
dclocalized b electrons, there is a loss of polaron energy 
proportional to this number, i. c.. 



The estimation of the kinetic energy entails probably 
the crudest approximation of this analysis. We assume 
that the the density of states is "flat" 



5(e) 



1 



6t< e<6t 



(17) 



Ep{nb, a, L) = EjTUbR^ 



(11) 



With this assumption we estimate the kinetic energy (for 
a spherical clump) to be 



Finally, there is the dclocalized kinetic energy of the elec- 
trons, which is expected to be of the form 



EKE{nb,R) 



An 



i?3(6i)nfc(l-nfc) 



(18) 



EKE{nb, R) = tlC{nb, R) 



(12) 



The total energy is given by 



where /C is a function of nb,R. It can be evaluated by 
finding the kinetic energy of bbR^ electrons in a clump of 
size R^. Note that we do not take E^e to go as 1/i?^ 
as might be expected from a free electron picture since 
the lowest energy that can be attained a tight binding 
scenario is bounded below by ~6t (for a cubic lattice). 
Putting everything together, we get 



E{nb,a,R) = — 



—Voil-nb-xfR^fia) + 
15 

EjTUbR^ - 6tnb{l - nb)R^) (19) 



Etotirib, a, R) = Essirib, a, R) + Ep{nb, R) 



We need to minimize this energy as a function of rib, 
a and i?; this will give us the clump size, spacing and 
charge distribution. 

In the remainder of this section R stands for a dimen- 
sionless number (clump size normalized by lattice param- 
eter a). Similarly charge densities are all dimensionless. 
We shall now proceed to estimate the different energy 
functions noted above. 

We estimate the electrostatic energy by assuming the 
box is "spherical" . In this case the total charge in the box 
is zero, and hence the electric field outside the sphere van- 
ishes, hence to a very crude approximation one box does 
not affect the energy of the neighbor. This is expected 
to be quite reasonable when the Ak ions are distributed 
homogeneously (as we have assumed in this calculation) , 
since the electrostatic interaction is expected to be well 
screened. With these assumptions, we find, 



The analysis of the minima of the above function shows 
that EjT < 6t is a necessary condition for the existence 
of a clump. If EjT > Qt, then the coefficients of R^ and 
i?^ are both positive (for any value of b) and hence R 
EKEi^b, R\15)Mii\ vanish for energy minimum. Inspecting ^ we see 
that this will result in a homogeneous distribution of the 
holes in the "bulk" equal to the cation charge density - 
the "bulk" becomes neutral. 

We now consider Ejt < Gt. Furthermore, 



1. The b electron density in the clump is determined 
by "chemical potential balance" . For the case of a 
flat band density of states, this will imply that the 
"chemical potential" is Ejt and 



-Ej 



-6t 



1 

12"r 



-de 



Ub 



rib 



1 - 



E 



JT 



6t 



(20) 



The factor a is determined from the condition that 
the hole density in the bulk vanishes (charge con- 
servation resulting in 



EEs{nb, a, R) 
where /(a 



(4^)^ 
45 

is the function 



Vo{l-nb-xfR''f{a) (14) 



Within this framework, 
obtain 



I - Ub 



we minimize (IT 



(21) 



r. t. R to 



/(«) 



1 



1 Aa^ - 5a3 -f 9a - 5) 



R' 



5(a- 



(a3 - 1)2 
l)\a + 2) 



9nfc(6t(l - m) - Ejt) 



(22) 



(a3 - 1) 



5(a- i; 



(15) 



= 6, 



The function / is such that /(I) = 1 and /(oo) 
a very slowly varying function of a and always of order 
unity. 



Note that the clump size varies as 1/ Vn) ■ 

We take Ejt ~ 2.5t, x « 0.3 to study the clump size 
as a function of Vq. The result is shown in FIG. [3l For 
these values, the centers of the clumps are spaced at a 
distance of about two and a half times the clump size 
(a = 1.33). A reasonable estimate of Vq in manganites 
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FIG. 3: Electron puddle size (normalized by lattice parameter 
a ) as a function of the Coulomb interaction parameter Vo 
(normalized by i) for Ejt ~ 2.5t, x = 0.3. 

is about O.li (Vb is likely to be between 0.01 and 0.1, see 
discussion below). In this range we see that the clump 
size is between 5 and 10 lattice spacings (taking a = 5 
A, we get i? = 2 - 5 nm), in surprisingly good agreement 
with available experiments. We must emphasize that we 
really do not have tight control of the constants (such 
as 9/47r appearing in the equation for i?). However, the 
main point that we learn is that the clump size is a few 
lattice spacings, much as what the Coulomb interaction 
is expected to do. We note here that there has been a 
recent calculation^^ of "phase separation" in doped man- 
ganites adding in Coulomb effects to the two fluid i — h 
model of Ref. [l^ along lines similar to those discussed 
above. However, that calculation does not uncover the 
dependence of the clump size on the long ranged Coulomb 
interaction as we have done here. 

We note again that the key assumption of the above 
analysis is the homogeneous distribution of the Ak ions. 
A more realistic solution that takes into account the in- 
homogeneous random distribution of the Ak ions requires 
a full scale numerical treatment of (U), and wc turn to 
this next. 



D. Approximate determination of the ground state 

The key approximation that permits a numerical de- 
termination of the ground state of the Hamiltonian ^ 
is a Hartree-like treatment of the Coulomb interactions 
between the h electrons. This allows us to treat all the 
electrostatic energy contributions in a classical fashion, 
i. e., replace the site charge operators qi by their expec- 
tation values qi in the ground state | ) : 

q, = l,\h%-h\h\). (23) 

In the absence of h states the problem reduces to the 
Coulomb glass problc m^'^i'^^ - a fully classical problem. 



The presence of the h states makes the system a new 
kind of quantum Coulomb glass - the h — b glass - with 
coexisting localized and delocalized states, in contrast to 
usual quantum Coulomb glasses^. Since the Ak ions are 
randomly distributed, and in addition the density of b 
electrons is very low, a Hartree approximation is likely 
to be reasonably accurate. Wc do indeed find (see be- 
low) that the interaction is well screened and the charge 
distribution is almost homogeneous at scales larger than 
a few lattice spacings. 

Our method of determination of the ground state 
is based on a generalization of the method previously 
used in finding ground states of the classical Coulomb 
glass^^i^. Here, the ground state is obtained by starting 
from a trial configuration (usually a random state), and 
"performing" transfers of electrons that lower energy un- 
til no transfer is possible that can lower the energy. In 
the classical Coulomb glass each transfer involves mov- 
ing an electron from its current position to a vacant site, 
creating a particle hole "excitation" . To follow a strat- 
egy similar to the one above in the present context, we 
investigate the energetics of possible transfers in the h — b 
glass. 

The following definitions are useful to understand ex- 
citations in the h — b glass. Each clump a has many 
'delocalized' one particle levels (obtained by diagonaliz- 
ing the kinetic energy on this clump); we label these one 
particle levels by r, and denote their kinetic energy by 
ea.r ■ The operator 

^L^E^r^l (24) 

i{a) 

creates a b electron in the r-th band state of the clump a, 
where ''^ is the associated single particle wave function 
spread out over the clump {i runs over all hole sites in 
the clump a as indicated by i{a)). Clearly, for a hole at 
site i 

g, = l-^|*rP^(M-ea.) (25) 

Here a is the clump to which the hole i belongs, and r 
runs over all the occupied band states in the clump a. 
The potential (pi at hole site i is defined as 

<f>^^VoY, — qJ. (26) 

Consider, for a given x, a configuration containing Nh 
holes and N}, band electrons (obviously, J^i 1i = — 
Ni, = xN), with Nc clumps. The extra energy in the 
system due to the addition of a hole at a hole- vacant site 
i without changing the clump structure is, 

E'l + + (27) 

The addition of a hole can, in general, change the clump 
structure since the new hole can modify existing clumps 
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or create new clumps. If the modified clumps contain 
band electrons, then their energies will change. Wc ignore 
these effects in writing ([27]) . A similar definition of ^, 
the energy increment for the addition of a band electron 
in the vacant band state a, r, is 

i 

+Vo ^m-^l'lK'^l', (28) 

where i,j run over all the hole sites in clump a. The 
definitions and E^ ^ serve as the equivalent of the 
"single particle levels" in the present context (Hartree 
approximation) . 

These definitions allow us to determine energies of "ex- 
citations" in the h — b glass corresponding to the cre- 
ation/annihilation of holes or 6/£-electrons. Strictly, the 
clump structure changes in an excitation; as mentioned 
earlier we ignore this and our excitations are "frozen- 
clump excitations" . 

Three types of excitations are possible in the h — b 
glass: 

1. h" excitations: Here a hole-occupied site j ob- 
tains an electron from a hole- vacant site (£ occupied 
site) i. The energy for this excitation is 



This type of excitations can affect the clump struc- 
ture. This arises from the fact that the clump struc- 
ture is determined by the position of the holes. 

2. "£ — 6" excitations: Here a hole-vacant site i (con- 
taining an £ electron) , donates an electron to a va- 
cant band state (a,r). The energy of this type of 
excitation is 




where j run over all the holes-sites in clump a. 
Again, this type of excitation can change the clump 
structure. It must be noted that the reverse of this 
process is also a possible excitation {a. b — £ exci- 
tation), in that a band electron annihilates a hole, 
the energy of which is negative of ([50)1 . Again, since 
both £ ~b and b — £ excitations change the position 
of the holes, they can affect the clump structure. 

3. "6 — 6" excitations: Here a 5-electron transfers from 
{a, r) b state to the vacant s) b state. The en- 
ergy of this excitation is 

^t.;a,. = Els - El, - E ^\^^'r\^r\"^ (31) 

t,3 



where i,j run over f3 and a (the last term requires 
obvious modifications if a = f3). Note that this 
type of excitation will not affect the clump struc- 
ture since it does not affect the position of the holes. 

For a given doping, the ground state is obtained in 
two stages. In the first stage, b states are not accounted 
for and the classical Coulomb glass ground state of the 
holes in presence of the electrostatic potential from the 
Ak ions is obtained. This is achieved by performing a 
scries oi £ — h excitations until a minimum energy con- 
figuration is achieved. In the second stage, starting from 
the classical Coulomb glass state, all possible excitations 
are performed iteratively. Each iteration consists of the 
following steps: 

1. Find best excitation: From the definition of the sin- 
gle particle levels E^, E^ and the excitation en- 
ergies defined above, the best possible excitation 
(the one that reduces the energy the most) is de- 
termined. 

2. Perform the excitation: If it is an £ — ft, excitation, 
then one hole is removed and one is added. If it is 
an £ — b or b — £ excitation, then either a hole is 
removed or added. If it is — 6 excitation, holes 
do not change. Thus, £ — h, and £ — b excitations 
change the clump structure. 

3. Update clump structure: If the clumps with occu- 
pied band states are not disturbed, then no update 
is required. If some of the clumps have been dis- 
turbed, then the intersection of the old clumps and 
new clumps is found, and electrons are distributed 
in the new clumps so that the charge distribution is 
as close as possible to that of the previous distribu- 
tion. After update, the number of band electrons 
must be same as before, unless isolated holes in the 
new clump structure are annihilated by band elec- 
trons. All energies are recalculated on update. 

Iterations are carried out until an energy minimum is 
achieved. It must be noted that the ground state that 
we obtain is based on single particle excitations. In gen- 
eral, the stability of the ground state must be checked 
for two (and multi) particle excitationsi^ This process is 
computationally intensive. The correctness of the ground 
state obtained in our case is affirmed by the energies of 
the highest occupied states of the £ electrons and the 6 
electrons - if the ground state calculation is correct, then 
the highest occupied states of both types of electrons 
will correspond to the chemical potential fi. We have 
never found serious violation of this criterion in several 
ten thousand simulations; we therefore believe that our 
scheme is robust enough to determine the ground state 
of the h — b glass. 

The electrostatic energy is calculated using the Ewald 
technique^^ using fast Fourier transform routines made 
possible by the use of periodic boundary conditions. The 
other computationally intensive step is the calculation of 
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FIG. 4: (color online) Density of states (DOS) of I polarons and 6 electrons, where their respective energies are defined by 
(|27|l and (|28|) . The occupied states are shaded. The chemical potential ^ and the center of the 6-band Eh are indicated. The 
difi'erent panels (a)-(d) show the variation of the densities of states as a function of Vo for Ejt = 1.0 and x = 0.3. The polaron 
density of states shows a Coulomb gap at the chemical potential for all values of Vb including Vo = 0.01 (graph (a), top- left). 



energies and wavef unctions of the h states in a clump. 
This is the most CPU intensive step which limits the size 
of the simulation cell. The size of the cell that can be 
calculated depends on the doping. For x ~ 0.2 we have 
calculated with cells as large as 20x20x20. 

Some other points regarding our simulations may be 
noted. First, the strong correlation U that induces many 
body effects is treated almost exactly (since U in the 
real system is very large, the U = oo limit is accurate); 
the 5-electron quantum dynamics is treated almost ex- 
actly (6-electrons do not hop to sites with t polarons). In 
our calculation, the long ranged Coulomb energy (within 
the Hartree approximation) is treated accurately using 
Ewald techniques (as mentioned above), and no further 
approximations are made in the calculation of the long 



ranged electrostatic energy. Further, we note that the 
Hamiltonian we study and the method we develop here 
to find the ground state is a novel generalization (both 
model and method) of the quantum Coulomb glass. The 
key point is that our model consists of two types of elec- 
tronic states, one localized polaronic (classical) and other 
delocalized band- like (quantum), and our treatment ac- 
counts for both of these on an equal footing. 



III. RESULTS 

This section contains the results of our study of the 
extended £b Hamiltonian ^ . The key parameters in the 
Hamiltonian are energies t, Ejt,Vo, and the doping x. 
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FIG. 5: (color online) Density of states (DOS) of £ polarons and 6 electrons, where their respective energies are defined by 
(|27|l and (|28p . The occupied states are shaded and the chemical potential is marked by fi. The center of the 6-band Et is 
also indicated. The difi'erent panels (a)-(d) show the variation of the densities of states as a function of x for Vo = 0.1 and 
EjT = 2.5. It is evident that the b bandwidth increases with increasing doping x. 



The hopping amplitude t is taken as the basic energy 
scale, and Ejx and Vq henceforth stand for the dimen- 
sionless values of the Jahn- Teller energy and the long 
range Coulomb interaction strength parameter normal- 
ized by t. Further all length scales are normalised by the 
lattice parameter a. We have checked for the size de- 
pendence of the results and found that results for cubes 
larger than 8x8x8 are essentially independent of the 
size. We checked densities of states, clump size distribu- 
tion etc., calculated in some 

cases up to cubes of size 20 x 20 x 20. We found that 
results followed essentially the same trends, independent 
of size. All the results shown here are averages from 
calculations obtained with one hundred random initial 
conditions using 10 x 10 x 10 cells unless stated otherwise. 



The density of states obtained from simulations is dis- 
cussed first. The energies of the "single particle states" 
are defined in (|27p for a hole (with a similar expression 
for an £ polaron) and (|28|) for the band electron. Figures 
mo and [5] show plots of densities of states for various val- 
ues of Vo, EjT and x, with the chemical potential fi and b 
band center Et indicated in each case. FIG. H] shows the 
variation of the density of states with the parameter Vq, 
FIG. [6] shows the variation in the density of states with 
the parameter Ejt (Vq and x fixed) and FIG. [5] shows 
the effect of doping on the density of sates for fixed values 
of Vq and Ejt- From a study of these figures we observe 
the following: The chemical potential is given by 

/i = -Ejt + Aix)Vo (32) 
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FIG. 6: (color online) Density of states of (DOS) £ polarons and 6 electrons, where their respective energies are defined by 
(|27|l and (|28p . The occupied states are shaded and the chemical potential is marked by fi. The center of the 6-band Et is 
also indicated. The different panels (a)~(d) show the variation of the densities of states as a function of Ejt for Vo — 0-1 and 
X = 0.1. The main effect of Ejt is to determine the chemical potential. Ejt does not affect width of the polaron density of 
states. The bandwidth of the b band is not strongly affected by Ejt, while the occupancy is strongly affected. 



where A{x) is a "very weak" function of x, with size 
of order unity; the physics behind this will be discussed 
later. 

We now discuss the £ polaron states:- As discussed in 
Sec. Ill Bl the polarons with long range Coulomb interac- 
tion form a Coulomb glass^^. A Coulomb glass possesses 
a "soft gap" at the chemical potential ^, in that the den- 
sity of states scales as 



Pi{e - ^) - (e - ^if 



(33) 



as is seen in FIG. |4l It is also evident that the width 
(energy spread) of the two lobes in the polaron density 
of states separated by the chemical potential increases 
with increasing Vq. Thus the polaron energies do not 



have a single value —Ejt, but have a distribution whose 
width scales with Vq. The energy spread arises out of 
fluctuations in the local electrostatic potential. This is 
the basic reason why though the polarons form a narrow 
band, there are no heavy fermion like specific heat ef- 
fects. The occupied polaron states (indicated by shaded 
portion below the chemical potential in FIG. [5]) per unit 
volume are proportional to (1 — x) for large values of 
Vo{> 1.0). For lower values of the Vq, the occupied spec- 
tral weight is less, owing to the fact that some electrons 
arc promoted to the b states. The effect of Ejt on the 
polaron density of states is shown in FIG. [H Since Ejt 
determines the chemical potential via ([5^ . it does not 
affect the width of the distribution, but only affects the 
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FIG. 7: (color online) Position correlation function for h — h and h — Ak. The figures in the left column are for random 
distribution of Ak ions, while those in the right column are for uniform distribution of Ak ions (uniform distribution means 
that the total charge of the Ak ions is distributed equally among the Ak sites) . The result is from simulations with a 20 x 20 x 20 
cube (for a single realization of the random distribution of the Ak ions in the left-side graphs). 



£ occupancy (amount of spectral weight in each lobe). 

Next, we turn to the b electron density of states. Sev- 
eral points may be noted. First, it is clear from FIG. [4] 
that for Vb < 1.0, the b band is centered at Et ~ B{x)Vo 
{B{x) is of order unity). The b band center is indepen- 
dent of EjT (see FIG. E}. Second, as is evident from 
FIG. [5l the bandwidth of the b states increases with in- 
creasing X (in fact, as ~ ^/x as we shall discuss below), 
and again essentially independent of Vq (FIG. [4]) as well 
as of EjT (FIG. [5]). Finally, the number of occupied b 
states (shown by shaded region in the figures), decreases 
with increasing Vq and Ejt- 

The physics underlying these observations may be un- 
derstood by studying the real space structure of the 
ground state, for example by studying the positional cor- 
relation of the holes present in the system (including the 
holes that appear due to promotion of the £ polarons to 
b electrons as well as the ones already present due to 
doping). FIG. [7] shows both the h-Ak and the h — h cor- 
relation functions, which measure, respectively, the prob- 
ability of finding a hole at a relative distance r from an 



Ak ion, or from another hole. (Note that the Ak ions are 
placed randomly in the A sites of the perovskite lattice.) 
If the system were purely random, the probability for any 
given value of r will be equal to the doping level x. When 
Vq is small (Vq = 0.01 in FIG. [T]), we see the hole-hole 
correlation function reaches a plateau at a distance r > 2 
with a value larger than x. This latter is because of the 
increased number of holes in the system due to £ polaron 
to b electron promotion. At larger Vq, the plateau in the 
hole-hole correlation function appears at x. The hole-Ak 
ion correlation function reaches a plateau of x, indepen- 
dent of Vo, as expected. The most important feature of 
the correlation function appears at r < 2. We see that the 
probability of finding a hole in the neighboring shells is 
non-zero and is highest at r = -s/S- Further, the probabil- 
ity of finding the Ak ions near the hole is also increased, 
though the increase is smaller than that of finding a hole. 
This suggests that there is a natural clustering tendency 
in the h — b glass where the holes tend to cluster around 
Ak ions. Clearly, this is due to their opposite charge 
causing gain in electrostatic energy. However, the com- 



13 



3 - 



2 



■ Simulation 



-DMFT(U = oc) 



0.1 



0.2 0.3 

X 



0.4 



0.5 



FIG. 8: The effective fialf bandwidth D^f / obtained from sim- 
ulations compared with the DMFT prediction ([2]). The half 
bandwidth in the simulation is calculated by using the clump 
structure obtained by minimizing the energy of the Coulomb 
glass. One hundred initial configuration are averaged over to 
obtain the simulation curve with the standard deviation bars 
indicated. 



pcting, repulsive electrostatic energy between the holes 
contributes to control this clustering tendency. Since the 
distance between the Ak ions and holes is smaller than 
that between two £ polarons, the clustering of holes near 
Ak ions will offset the energy penalty of £ polarons com- 
ing together - the short (Angstrom sized) clustering scale 
is result of this competition. Strong screening effects are 
evident in the correlation functions of FIG. [7]- the cor- 
relation functions reaches a plateau at about 2-3 lattice 
spacings. 

The observations above allow us to obtain the chemical 
potential of the system. The hole clustering caused by 
the hole-Ak ion interaction causes an additional effective 
electrostatic potential at an £ electron site which can be 
written as A{x)Vo where A{x) is a (here undetermined) 
doping dependent factor of order unity. The quantity 
A{x)Vo can be interpreted as the average electrostatic 
potential at £ polaron sites. The actual energies at the 
sites fluctuates about the mean total energies of the £ 
polarons given by —Ejt + A{x)Vo- From the Coulomb 
glass problem it is known that the chemical potential will 
be equal to the average energy^^i^ and ([32]) follows. 

We now discuss the energetics and the nature of the b 
states, since the transport properties depend crucially on 
the nature of the b states and whether they are occupied. 
In case they are not occupied or if all the occupied b 
states are localized, the system is an insulator; on the 
other hand if extended b states are occupied one has a 
metal. We summarize our results in this regard, and 
compare them where ever possible with those obtained 
in single site DMFT-iiii^iiS 

The hole site has an average electrostatic potential en- 
ergy B(x)Vo where B(x) « 1; this follows from argu- 
ments very similar to those used earlier for the electro- 



static (Hartree) energy of the £-polaron sites. Since the b 
electrons occupy clumps consisting of hole sites, they too 
sense this average repulsive potential B{x)Vo and thus 
the band center Eb shifts to this value (as is seen in 
the simulation, see FIGSHl El El). The simulations in- 
dicate that the effective bandwidth -De// of the b states 
depends only on the doping x. To understand this re- 
sult, we calculated D^ff as a function of x using clumps 
obtained from a Coulomb glass calculation of the holes 
including the random distribution of the Ak ions (thus 
the calculation was performed with only £ — h excita- 
tions). The result of this calculation is shown in FIG. [5] 
along the the DMFT result for the bandwidth. It is evi- 
dent that the agreement is excellent for a wide range of x. 
Note that the half bandwidth is provided by the clump 
structure determined by long ranged Coulomb interac- 
tions and the random distribution of the Ak ions which 
induce the clustering tendency of holes around them (as 
discussed above). Thus the result for the effective band- 
width (21) seems to arise out of two key factors: (a) The 
large U limit which disallows simultaneous b and t occu- 
pancy, and (b) The clumps of holes induced by the ran- 
dom distribution of the Ak ions. This is interesting since 
it implies that the single site DMFT which neglects long 
range Coulomb interactions and consequent nanoscale in- 
homogeneities or clumping, and replaces it with a self 
consistent annealed random medium in which all sites 
are equivalent, is accurate for such purposes. We also 
note here that in the DMFT, the £ polaron is a single 
level while in the simulations the occupied levels, though 
localized, have a distribution of energies. 

We now discuss the doping/hole density Xci at which 
b states are occupied. In the DMFT, for [/ = oo and at 
T = 0, the critical concentration Xc at which the lowest 

energy b state is occupied is Xd = \ — — ) where Dq is 

the bare half bandwidthi^. From our simulations, the £ 
chemical potential (i. e., the energy of the last occupied 
£ state) is —Ejt+A{x)V[). At the onset of the occupation 
of the b states, equilibrium demands that the bottom of 
the b band has to just equal the chemical potential ji. 
Noting that the effective half bandwidth of the b band is 
Df)^/x and the band center E^ = B{x)Vq^ one obtains on 
equating the b band bottom to the chemical potential, a 
result for the critical doping 

Vi^T = ^ + {A{x,) - B{x.)) ^. (34) 

Since A{x) « B{x) « 1, and Vo/Da < 0.05 (typically) 
Xci is expected to be very close to the DMFT value. In 
FIG. El we notice that as expected, for small/realistic 
values of Vb, there is good agreement with DMFT and the 
simulation results. For larger values of Vb (of the order of 
Dq), b state occupancy needs much larger hole density x 
than implied in the DMFT. Indeed, there is critical value 
of Vb(Vj3'^ w 1.2) above which there is no band occupancy 
for Ejt > 1; the system becomes a Coulomb glass of £ 
polarons. 
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FIG. 9: (color online) Critical doping levels Xd and Xc2 obtained from simulations. The lightest region in the contour plot 
contains no b electrons, the intermediate shade has "trapped" 6 states occupied, and the darkest regions corresponds to b states 
that percolate through the simulation box. The solid line corresponds to the analytical DMFT result for Xci- 



Furthermore, unlike in the DMFT, in the present con- 
text the fact that b states are occupied does not in itself 
make the system a metal. Even when x > Xci, the hole 
sites where b states exist form compact clumps in our sim- 
ulations. The b levels within the bigger clumps which are 
below the chemical potential and are occupied are hence 
still localized. FIG. [11] shows the 'clumps' with both 
the occupied and unoccupied b state regions indicated. 
We notice that the clumps, and hence the b puddles are 
generally isolated. Hence the system is still an insula- 
tor. The system becomes metallic for x > Xc2, when the 
6 puddles percolate, which requires as a necessary condi- 
tion that the clumps with the b states occupied percolate. 
Xc2 can be estimated by calculating the inverse participa- 
tion ratio^^ and by checking for the geometric percolation 
of the clump. The differences among Xci (DMFT), Xd 
(simulations), Xc2 (simulations) or in the nature of the 
b states and their occupancy, are due to the fact that 
the single site DMFT of the ib model (P) and the nu- 
merical simulation of the extended £b model (|3]) are very 



different approximations. In the former, the static £ po- 
larons are repulsive (potential U) scattering centers for 
b electrons, present randomly and independently at each 
site with a self consistently determined probability close 
to (1 — x). The scattering is treated in CPA (coherent 
potential approximation), so that there is no Anderson 
localization of the b electron states; there are diffusive 
and extended. By contrast, in the simulation, the b elec- 
trons are completely excluded from £ polaron sites, and 
can hop freely from site to site only among contiguous 
hole sites (clumps). So unless the clumps are are con- 
nected percolatively, the system is an insulator. Further- 
more, even when the clumps percolate, if the £ polaron 
distribution and the b electron propagation in such an £ 
polaron medium could be treated realistically, one will 
find that in contrast to single site DMFT, a fraction of 
the occupied b states are Anderson localized due to the 
£ polaron disorder and the repulsion Ununu- Hence Xc2 
is larger than the value of x at which the clumps just 
percolate, as the highest occupied b state could still be 
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FIG. 10: (color online) DC conductivity a{ui ~ 0) as a function of system parameters. For long long range Coulomb parameter 
Vb greater than 0.4, the conductivity is zero for the indicated range of Ejt and x. The darkest shade of the contour corresponds 
to a dimensionless conductivity a > 10, the next darkest 5 > cr > 10, the lightest > a > 5. 



Anderson localized within the percolating clump; hence 
signifies the value of x at which the mobility edge in 
the b band crosses the chemical potential. 

We have confirmed this by calculating the contribution 



of the b states to the DC conductivity using the standard 
Kubo formalism^! - the results, in the form of contours of 
constant dc conductivity in the Ejt, x plane, are shown 
in FIG. [TOl Again, it is clear that the metal-insulator 
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(d) 



FIG. 11: (color online) Real space structure of the electronic state. The darkest regions (magenta in online version) denote 
hole clumps with occupied 6 electrons, the lightest (white) denote hole clumps with no b electrons, the second lighter shade 
(cyan) denote singleton holes, and the second darkest shade (light blue) represents regions with £ polarons. The simulations 
for each Vb are for the same realization of the random distribution of Ak ions. The cell size is 16 x 16 x 16. 



boundary (the Xc2) curve is independent of the long range 
Coulomb parameter Vq when Vq is small. However, the 
DC conductivity completely vanishes for larger values of 
Vq greater than V^'" of about 0.6 and the system is an in- 
sulator (due to both strong local disorder potential from 



the Ak ions and Coulomb repulsion). Not unexpectedly, 
this critical value Vq^ is much smaller than the value Vq 
that forbids b state occupancy. 

Finally, we discuss the clump size distribution. The 
size of a clump is defined in terms of a radius of gyra- 
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FIG. 12: (color online) Size distribution of clumps with occupied b electrons. Clump concentration is number of clumps per 
lattice site (simulated with a 10 x 10 x 10 box). Size of the clump R is the "radius of gyration" of the clump (see text), ordinate 
is inverse clump size. 1/ii = corresponds to a percolating clump, while 1/R = 2 corresponds to the smallest clump with two 
sites. The plots show the effect of the Coulomb interaction Vb on the distribution for x = 0.3 and Ejt = 2.5. For Vb = 2.00 
there are no clumps with occupied 6 electrons. 



tion R calculated as i? = 




r)2 where r 



is the center of mass of the clump and i runs over all the 
holes in the clump. We investigate the distribution of 
the quantity (l/R) as a function of the system parame- 
ters for clumps with b states occupied. For a percolating 
clump the quantity 1/R and for the smallest clump 
with two holes 1/R = 2.0. Figures [T2l [T3l and [Til show 
the concentration (number per lattice site) of the occu- 
pied clumps as a function of their inverse size. For a given 
doping X and Ejt, there are many more larger clumps 
for smaller Vo (FIG. \1^. At high enough Vb (Vb = 2.0) 
there are no b electron puddles. The effect of doping x 
is also as expected. At small doping, there are clumps 
of various sizes, while at larger doping, there is only one 



percolating clump (see FIG. [T3|) . The effect of Ejt is 
more interesting (FIG. [Ti|) . For Ejt — 1.50 there is a 
percolating clump together with isolated small puddles. 
On increase of Ejt {Ejt = 2.0), the isolated electron 
puddles are larger in size, and on further increase of Ejt 
{Ejt = 3.0), there is essentially a single large percolating 
puddle. This is due to the fact that occupancy of a clump 
requires that the gain in kinetic energy offsets the loss 
in the Jahn- Teller energy Ejt at the least, thus larger 
Ejt requires larger clumps to be occupied for sufficient 
gain in kinetic energy. FIG. [15] shows the dependence 
of the average inverse size {1/R°'^)ave on the parameters; 
the results are as expected, in that the average size of 
clumps is about a few lattice spacings for small doping 
{x < 0.25), followed by a regime of percolating clumps. A 
result of particular interest is the dependence of the aver- 
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FIG. 13: (color online) Size distribution of clumps with occupied b electrons. Clump concentration is number of clumps per 
lattice site (simulated with a 10 x 10 x 10 box). Size of the clump R is the "radius of gyration" of the clump (see text), ordinate 
is inverse clump size. 1/i? = corresponds to a percolating clump, while 1/R = 2 corresponds to the smallest clump with two 
sites. The plots show the effect of the doping x on the distribution for Vb = 0.1 and Ejt = 1.0. For x > 0.4 there is only one 
percolating clump. 



age clump size on the long range Coulomb parameter Vq. 
This result is plotted in FIG. [1^1 we see that the clump 
size is essentially insensitive to Vq until Vq approaches 
t. This may be contrasted with the simple analytical re- 
sult obtained earlier (see FIG. [3]). It is evident that the 
the real system, any finite Vq will prevent "phase separa- 
tion" , but the final details of the ground state is strongly 
determined by the random distribution of the Ak ions, 
which also determines the size scale of the clumps (and 
hence their insensitivity to Vq for small values of Vb). We 
have investigated the size dependence of the clump size 
on the size of the simulation cube and found that there 
is no significant size dependence. 

Our study therefore indicates that for small doping 
(x < 0.25), electronic inhomogeneities in the the form of 
localized b electron "puddles" exist in a background of 



polarons (see FIG. [TT|) while for larger doping, the sys- 
tem consists of intermingled states of percolating, delo- 
calized b electron puddles coexisting with regions of local- 
ized £ polarons. The scale of the inhomogeneities is that 
of the local Ak ion disorder and not much larger. The 
study therefore confirms that long range electrostatic in- 
teractions gives rise to only nanometer (or lattice scale) 
electronic inhomogeneities as was anticipated in earlier 
literaturei^ The key point is that this nanometric scale 
arises not out of phase competition, but due to strong 
correlation effects of two types of states that appear at 
the atomic scale, whose spatial arrangement is then con- 
trolled by the long range Coulomb interaction and the 
random distribution of the Ak ion disorder. 
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FIG. 14: (color online) Size distribution of clumps with occupied b electrons. Clump concentration is number of clumps per 
lattice site (simulated with a 10 x 10 x 10 box). Size of the clump R is the "radius of gyration" of the clump (see text), ordinate 
is inverse clump size. 1/ii = corresponds to a percolating clump, while 1/R = 2 corresponds to the smallest clump with two 
sites. The plots show the effect of Ejt on the distribution for Vb =0.1 and x = 0.3. For Ejt = 3.5 there are no occupied 
clumps. 



IV. DISCUSSION 

Our study of the extended ib model has established 
that long ranged Coulomb interaction gives rise to 
nanoscale electronic inhomogeneities. A key point is that 
even a very small Vq completely eliminates the macro- 
scopic phase separation tendency of the Falicov-Kimball 
like £b model (compare FIG. [1] with FIG. [11]). By using 
a dielectric constant of about 20, we estimate an upper 
bound for Vq in manganites to be about 0.02 eV and thus 
well within the small Vq regime; for t sa 0.2 eV, the di- 
mensionless Vq is about 0.1. Furthermore, for a given 
set of energy parameters (Vb and Ejt), we have shown 
the existence of two important doping thresholds: Xd 
which corresponding to the occupancy of the b states, 
and Xc2 > Xci at which there is an insulator to metal 



transition. Our results for Xd are close those of the ear- 
lier DMFT prediction ; however, the existence of Xc2 is 
a new aspect of this work (in DMFT Xc2 — Xci)- The 
key physics behind the agreement of Xci is the large U 
condition (no simultaneous occupancy of £ polaron and b 
electron on one site), and the random distribution of the 
Ak ions. As we have shown, the random distribution of 
Ak ions produces local clustering of holes around them 
which eventually causes electron puddles and subsequent 
percolating clumps at higher doping. The main conse- 
quence is that the electronic inhomogeneities are all at 
the nanoscale and the material appears homogeneous at 
microscalc. 

How important is the random distribution of Ak ions? 
This question was answered by distributing the charge of 
the Ak ions equally among all the Ak sites. The physics of 
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FIG. 15: (color online) Average inverse size {1/ R°)a of the clumps with at least one b state occupancy. Clump size R is defined 
as the radius of gyration (see text) of the clump. For small Vb, the average inverse size at a; = 0.2 and Ejt = 1.0 is between 
0.44 and 0.88 - thus corresponds to spherical puddle of containing 6 to 10 lattice sites. With increase in Vo, the sizes decreases 
and eventually there are no occupied clumps at Vb = 2.0 (not shown). At larger x, there is one percolating clump for Vb < 1.0. 
At larger Ejt, again, there are no occupied clumps. 
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FIG. 16: Clump size as a function of the Coulomb interaction 
parameter Vo (normalized by t) for Ejt = 2.5t, a; = 0.3 
obtained from simulations. 



the problem is changed completely by this step. First, the 
clustering tendency of the holes is completely suppressed. 
This is evident from FIG. [7] (right column) where it is 
seen that the probability of finding a hole neighboring a 
hole is nearly reduced to zero. In such a case, for a given 
set of energetic parameters, much larger doping levels 
would be required for b state occupancy, and for the insu- 
lator to metal transition. Indeed, these observations are 
borne out in the full scale simulations, see FIG.[T71 Thus 
the physically more realistic random distribution of Ak 
ions is key to the agreement with the DMFT calculation. 
In other words, the earlier DMFT calculation, although 
it did not take into account the long range Coulomb in- 
teractions, premised a homogeneous state (on a "macro- 
scopic" scale), thereby effectively incorporating the key 
effect of the long ranged Coulomb interaction. 

It may be argued that effects of the random distri- 
bution of Ak ions (that help the clustering of holes) is 
equivalent to a local disorder potential. The results of 
the £b Hamiltonian ^ with an additional disorder po- 
tential (but no long range Coulomb interactions). 



(35) 



where Wj are distributed uniformly between — W and W ^ 
should be roughly similar to those of the extended Ih 
Hamiltonian ([T]). The premise is that since I electrons 
will occupy sites of low Wi , holes are likely to clump with 
a probability proportional to where n is the size of the 
clump. This is similar to the probability of finding clus- 
tering together of Ak ions. Thus the clumps that appear 
in the disorder only model may be expected to have sim- 
ilar nature as that of the case with random distribution 
of Ak ions. The results of the disorder-only Hamiltonian 
are shown in FIG.[TS1 Clearly, the result of Xci is close to 
that of the DMFT and the extended £6 results. Also, the 
clump sizes etc. are also quite similar as expected. There 



is, however, a crucial difference - the density of states of 
the £ electrons does not have a Coulomb type gap in the 
purely disordered model just discussed. In the real sys- 
tem both disorder and long ranged Coulomb interactions 
are present, and the extended (. — b model includes both 
in a realistic way, in contrast to the disorder only model. 

We now discuss further predictions and inferences from 
the model. Our model is particularly suited to study the 
low bandwidth manganites such as Pr based compounds 
which have a large fcrro-insulating regime as a function 
of doping. Most of the inferences made hinge on the 
coexistence of localized and delocalized states - the key 
physical idea of the model. It may be inferred from the 
extended £6 model that the low temperature conductiv- 
ity in the low temperature regime for doping x < Xci will 
be governed by the polaron Coulomb glass. Thus, the 
low temperature conductivity is expected to be that pre- 



dicted by for the Coulomb glas; 
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1. e. 



a(r) 



-a/Vr 



where a is a constant. Further, above Xc2 > x > Xci, sev- 
eral excitations contribute to the conductivity. These will 
include polaron hopping, variable range hopping of the 
6-state electrons from one puddle to the other, even two 
step processes where a 6-state in an intermediate for an £- 
polaron to hop from one site to another. Transport in this 
regime is therefore expected to be involved and further 
investigation is necessary to uncover the temperature de- 
pendence. We note that transport measurements^^ on 
doped manganites do show the features that we deduce 
from our model. For x > Xc2, the low temperature con- 
ductivity will be metallic with large residual resistivity. 
An important contribution to this resistivity arises from 
the Couloumb potentials of the Ak ions as has been noted 
earlier—. At higher temperatures, scattering from the 
thermally disordered t2g spins will cause a decrease in the 
effective bandwidth of the h states and can lead to the 
opening up of a gap between the £ states and the h states 
above the ferromagnetic transition temperature. In this 
high temperature regime the conductivity will have con- 
tributions both from polaron hopping and thermal exci- 
tations of the £ polarons to the h states similar to that of 
a semiconductor. 

The present model does not explicitly include the t2g 
core spins as degrees of freedom. A simple minded ap- 
proach of including these degrees of freedom will lead 
to intractable computational complexities. Novel ap- 
proaches to treat both Coulomb interactions, h state 
quantum effects, and core spins need to be developed 
to attain a final understanding of the manganite puzzle. 

We conclude the paper with a discussion of the im- 
portant issue of "phase separation" and electronic in- 
homogeneities. As noted in the introductory section 
several groups have implied that the presence of elec- 
tronic inhomogeneities is essential for the occurrence of 
the colossal magnetoresistance. Our work suggests that 
inhomogeneities are only present at the nanometer scale 
and a theory that averages over these without explicit 
treatment of these inhomogeneities can and does repro- 
duce colossal magnetoresistanceJii^ii^ Indeed, very re- 
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FIG. 17: (color online) Critical doping levels Xci and Xc2 obtained from simulations with uniform distribution of Ak ions. The 
contour scheme is same as that of FIG. [§] These results are for one initial realization of the holes (no average over initial 
conditions) . 




FIG. 18: Critical doping levels Xd and Xc2 obtained from simulations with on-site disorder which varies from —W to +W . The 
contour scheme is same as that of FIG. [9] These results are for one realization of the disorder. 
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FIG. 19: Schematic of a suggested experiment to study the 
influence of long range strain effects on mesoscale inhomo- 
geneities in manganites. 



cent experimental worki^ shows that materials with no 
detectable "phase separation" show colossal magnetore- 
sistance. Further, we find that nanoscale inhomogeneities 
are a direct result of long range Coulomb interaction frus- 
trating "phase separation" induced by strong correlation 
— nanoscale electronic inhomogeneities are not a result 
of phase competition in our model. Although similar 
mechanisms based on long range coulomb interactions 



has been discussed earlieriSiiSiiiidSiiM, we believe that this 
is the first detailed quantitative treatment of a realistic 
model any correlated oxide. 

The large length scale (micrometer sized) inhomo- 
geneities seen in experiments remain to be explained in 
the present framework. To the best of our knowledge, 
the coexistence of metallic and insulating regions have 
all been seen only in surface probe measurements or mea- 
surements with thin films (for electron microscopy). Cre- 
ation of a surface introduces defects such as cracks and 
steps all of which have long ranged elastic fields. As is 
well known, the phase of manganites is strongly influ- 
enced by pressure (stresses) 4^ Thus the large scale in- 
homogeneities are likely to be a result of pre-existing 
strain sources as is indicated by recent photo-emission 
experiments^. It is possible to explicitly test this in an 
experiment (see FIG. [T^ with a pre-cracked mangan- 
ite sample - on loading the cracked sample the motion 
of the metal-insulator boundary is expected to be ob- 
served. Similar experiments^ (without cracks etc.) do 
indeed suggest the strong effects of strains. 

A second possibility for the existence of micron-scale 
clusters could be due to "kinetic arrest" as seen recently 
in some rare-earth compounds^. Thus the patches that 
appear could arise to an "incomplete phase transition" , 
and likely to show "glassy" behavior. Indeed many man- 
ganites are known to show glassy behavior—, and this is 
possibly another important line of further investigation. 
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